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Abstract 

The algebraic structure underlying the totally asymmetric ex- 
clusion process is studied by using the Bethe Ansatz technique. 
From the properties of the algebra generated by the local jump op- 
erators, we explicitly construct the hierarchy of operators (called 
generalized hamiltonians) that commute with the Markov oper- 
ator. The transfer matrix, which is the generating function of 
these operators, is shown to represent a discrete Markov process 
with long-range jumps. We give a general combinatorial formula 
for the connected hamiltonians obtained by taking the logarithm 
of the transfer matrix. This formula is proved using a symbolic 
calculation program for the first ten connected operators. 
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1 Introduction 

The Asymmetric Simple Exclusion Process (ASEP) is a driven lattice gas 
of particles that hop on a lattice and interact through hard-core exclusion. 
Originally, the ASEP was proposed as a minimal model in one-dimensional 
transport phenomena with geometric constraints, such as hopping conductiv- 
ity, motion of RNA templates and traffic flow. The exclusion process displays 
a rich phenomenological behaviour and its relative simplicity has allowed to 
derive many exact results in one dimension. For these reasons, the ASEP has 
become one of the major models in the field of interacting particle systems 
both in the mathematical and the physical literature and plays the role of 
a paradigm in non-equilibrium statistical mechanics (for reviews, see e.g., 
Spohn 1991, Derrida 1998, Schiitz 2001). 

It has been shown that the evolution operator (or Markov matrix) of 
the exclusion process can be mapped into a non-hermitian Heisenberg spin 
chain of the XXZ type (Gwa and Spohn, 1992; Essler and Rittenberg 1996). 
This mapping allows the use of the techniques of integrable systems such as 
the coordinate Bethe Ansatz (for a review see, e.g., Golinelli and Mallick, 
2006b). Spectral information about the evolution operator (Dhar 1987; Gwa 
and Spohn 1992; Schiitz 1993; Kim 1995; Golinelli and Mallick 2005) and 
large deviation functions (Derrida and Lebowitz 1998) can be derived with 
the help of coordinate Bethe Ansatz. Besides, using the more elaborate 
algebraic Bethe Ansatz technique, the eigenstates of the Markov matrix can 
be represented as Matrix Product states over finite dimensional quadratic 
algebra (Golinelli and Mallick, 2006a). The algebraic Bethe Ansatz also 
plays a fundamental role in the derivation of the Bethe equations for ASEP 
with open boundaries (de Gier and Essler, 2005, 2006). 

The aim of the present work is to explore the algebraic properties of the 
totally asymmetric exclusion process (TASEP) that stem from the algebra 
generated by the local jump operators that build the Markov matrix. The al- 
gebraic Bethe Ansatz technique allows to construct a hierarchy of generalized 
hamiltonians that contain the Markov matrix and commute with each other. 
The generating operator for this family, called the transfer matrix, defines 
therefore a commuting family of operators that can be simultaneously diago- 
nalized. We derive, using the local jump operators algebra, explicit formulae 
for the transfer matrix and the generalized hamiltonians and characterize 
their action on the configuration space. These generalized hamiltonians are 
non-local because they act on non-connected bonds of the lattice. However, 
connected operators are generated by taking the logarithm of the transfer 
matrix. We study these connected operators and give an explicit formula for 
them. 
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The outline of this work is as follows : in Section 2, we describe the 
basic algebraic properties of the totally asymmetric exclusion process and 
define the associated algebra. In Section 3, we give explicit formulae for the 
transfer matrix and for the generalized hamiltonians that commute with the 
Markov matrix M. In particular, we show that the transfer matrix can be 
interpreted as a discrete time Markov process and we describe the non-local 
actions of the generalized hamiltonians. In Section 4, we study the connected 
operators obtained by taking the logarithm of the transfer matrix and propose 
a conjectured general formula for these local operators. The actions of these 
operators are described explicitly. Some mathematical proofs are given in 
the appendices. 

2 Algebraic properties of the TASEP 
2.1 Definition of the model 

The simple exclusion process is a continuous-time Markov process (i.e., with- 
out memory effects) in which indistinguishable particles hop from one site 
to another on a discrete lattice and obey the exclusion rule which forbids 
to have more than one particle per site. In this work, we shall study the 
case of particles hopping on a periodic 1-d ring (see figure [I]) with L sites 
labelled i — 1, . . . , L (sites i and i + L are identical due to periodic boundary 
conditions). The particles move according to the following dynamics: during 
the time interval [t, t + dt], a particle on a site % jumps with probability dt to 
the neighboring site % + 1, if this site is empty. This model is called 'totally 
asymmetric' because the particles can jump only in one direction. The exclu- 
sion rule forbids particles to overtake each other and their ordering remains 
unchanged. Moreover, as the system is closed, the number of particles is 
constant. 

The state of a site i is encoded in a boolean variable r», where n — 1 if 
i is occupied and = otherwise. The two-dimensional state space of the 
site i is noted Vi (we have Vi = C 2 ) and its basis is given by (|lj), |0j)). A 
configuration C of the system of L sites is written as 



The state space Tii of the ring is therefore a 2 L dimensional vector space 
given by 



Due to the conservation of the number n of particles, Hl splits into invariant 
spaces Ti^ 1 of dimension L\/[n\(L — n)\], characterized by X^iLi r « = n - 





H L = V x <g> V 2 ® . . . <g> V L . 



(2) 
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dt 




Figure 1: The totally asymmetric exclusion process on a ring. Sites are 
labelled from 1 to L; a particle jumps with probability dt to the neighbouring 
forward site if this site is vacant. 

The probability distribution of the system at time t can be represented as 
a vector ip(t) G Hl, where the component ipcif) is the probability of being 
in the configuration C at time t. The vector ip(t) evolves according to the 
master equation 



where M is a 2 L x 2 L Markov matrix acting on H L . For C ^ C\ M(C, C) 
is the transition rate from configuration C to configuration C: it is equal to 
1 if C is obtained from C by an allowed jump of a particle, and otherwise. 
The diagonal elements are negative and — M(C, C) is the exit rate from C, 
i.e., the number of allowed jumps from C. The sums over columns of M 
vanish, J2c M(C, C) = 0. This property ensures probability conservation : 



In the case of the TASEP on a periodic ring, sums over rows of M also 
vanish. This implies that the stationary probability, obtained for t — > oo, is 
uniform over each subspace . 

2.2 The algebra of jump matrices 

The Markov matrix can be written as 




L 




(4) 



i=i 



where the local jump operator Mi represents the contribution to the dynamics 
of jumps from the site % to % + 1. Thus, the action of the 2 L x 2 L operator Mj 
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affects only the sites % and i + 1 and is non-zero only if Tj = 1 and r i+1 = : 

M l \r 1 ...l,0...r L ) = |r 1 ...0,l...r L )-|r 1 ...l,0...r L ) , (5) 
M i \T 1 ...T i ,T i+1 ...T L ) = 0, if Ti = or r i+ i = f . (6) 

The operator M L corresponds to jumps from site L to 1 (the site L + 1 is 
identical to the site 1 because of the periodic boundary conditions). 

Using this definition of the local jump operators, it can be verified that 
the Mi satisfy the following relations : 

Ml = -M^ (7) 

Mi M i+1 Mi = M l+1 Mi M l+1 = 0, (8) 

[M,M,] = if > 1. (9) 

As we consider a periodic system, we use the convention Ml+i = Mi 
in the above relations. We emphasize that [M L ,M]\ ^ and M X M L M X = 
Ml Mi Ml = 0. 

The algebra generated by the Mj operators will be called here the TASEP 
algebra. We remark that the Mj operators can be obtained as a limit of the 
Temperley-Lieb algebra generators. We shall call, by definition, a word, 
any product of the Mj's; any element of the algebra can be written as a 
linear combination of words. The length of a word is the minimal number of 
operators Mj required to write it. 

Each word acts on the configuration space Hl and can be described as a 
series of jumps. For example, the word MiM 2 describes a jump of a particle 
from site 2 to 3, followed by a jump of another particle from site 1 to the site 2; 
the action of M X M 2 on a configuration vanishes unless t± — 1, r 2 = 1, r 3 = 
and we have 

M 1 M 2 |l,l,0,r 4 ...r L ) = |0, 1, 1, r 4 . . . r L ) - |1, 0, 1, r 4 . . . r L ) . (10) 

Similarly, the word M 2 M 1 represents a jump of a particle from site 1 to 2 
followed by a jump of the same particle from 2 to 3 : 

M 2 M 1 |l,0,0,r 4 ...r L ) = |0, 0, 1, r 4 . . . r L ) - |0, 1, 0, r 4 . . . r L ) . (11) 

Clearly, Mi and M 2 do not commute because the jumps on two adjacent sites 
are not independent. 

2.3 Ring-ordered product of jump matrices 

We define here the ring-ordered product of jump matrices which will be used 
in the following sections. 
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The ring-ordered product O () acts on words of the type 

W = M h M i2 ...M ik with 1 < i x < i 2 < . . . < i k < L , (12) 

by changing the positions of matrices that appear in W according to the 
following rules : 

(i) If ii > 1 or i k < L, we define O (W) = W. The word W is well- 
ordered. 

(ii) If ii — 1 and i k = L, we first write W as a product of two blocks, 
W = AB, such that B = M b M b+1 . . . M L is the maximal block of matrices 

with consecutive indices that contains M L) and A = \I\ \I; U, with 

i a < b — 1, contains the remaining terms. We then define 

O (W) = O (AB) = BA = M b M b+l . . . \I,.UM,, ■ ■ ■ M ia . (13) 

(iii) The previous definition makes sense only for k < L. Indeed, when 
k = L, we have W = M\M.2 ■ ■ ■ Ml and it is not possible to split W in two 
different blocks A and B. For this special case, we define 

<D(M X M 2 ...M L ) = |1,1,...,1)(1,1,...,1|, (14) 

which is the projector on the 'full' configuration with all sites occupied. 

The ring-ordering O () is extended by linearity to the vector space spa- 
nned by words of the type described above. 

Let us give some examples. For k — or 1, the ring-ordered product acts 
trivially: 0(1) = 1 and O (MA = Mj. For k = 2 , we have O (MjM,-) = 
MiMj when 1 < % < j < L; however, for the special case when % — 1 and 
j = L,0(M 1 M L ) = M L M 1 . 

The ring-ordered product embodies the periodic boundary conditions. On 
a ring, the natural order between integers is not valid. Indeed, Mi and Mi 
act as neighbouring bonds and site L should be viewed as being 'behind' 
site 1, just as site 1 is behind site 2. The ring-order product restores the 
correct order on a ring and allows to construct operators that are translation 
invariant. For example, for L — 3, the operator U = M l M 2 + M 2 M 3 + M l M 3 
is not well-ordered and does not commute with translations. But, O (U) = 
M\M 2 + M 2 M 3 + M3M1 is well-ordered and does commute with translations. 
Finally, we remark that when a ring-ordered product acts on a configuration, 
each particle advances by at most one lattice unit : indeed, because terms 
such as Mj +1 Mj do not appear in a ring-ordered product, no particle can 
perform multiple jumps. 
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3 Transfer matrix and generalized hamiltoni- 
ans 

The algebraic Bethe Ansatz is a method for diagonalizing the hamiltonian of 
integrable models (for a review, see, e.g., Korepin et al. 1993; for a pedagogi- 
cal introduction, see, e.g., Nepomechie 1999). This technique can be applied 
to the Markov matrix M of the TASEP (Golinelli and Mallick, 2006b). The 
key step is to construct a family of transfer matrices, t(\), which act on the 
configuration space TCl- For any value A and v of the spectral parameter, we 
have 



Thus, the operators t(X) form a one-parameter family of commuting opera- 
tors which depend on a real number A, called the spectral parameter. This 
family contains the Markov matrix M as will be shown below. Therefore, all 
the t(A)'s share with M a common eigenvector basis independent of A. For 
the TASEP, these eigenvectors are determinants of matrices involving the 
roots of the Bethe equations and the corresponding eigenvalues are functions 
of A (see e.g., Golinelli and Mallick 2006b for an explicit formula). 

The transfer matrix t(X) is a polynomial in A of degree L : we can thus 
define Hi,H 2 ,... H L as follows 



The iffc operators are 2 L x 2 L matrices acting on the configuration space H,^. 
The f/fc's will be called 'generalized hamiltonians' by analogy with quantum 
spin systems (Arnaudon et al. 2005). As the H^s are derivatives of t(X), 
they also commute with each other : 



for all j, k and A. More generally, any operator generated from t(X), or equiv- 
alently from the generalized hamiltonians, belongs to the same commuting 
family. 

The above considerations are familiar in the framework of algebraic Bethe 
Ansatz. In Appendix |A], we explain how the transfer matrix can be con- 
structed using this method. 

3.1 Expressions of the H^s 



[t(A),f(z/)]=0. 



(15) 




(16) 



[H k ,t(X)] = 0, [Hj,H h \={] 



(17) 



In this section we describe our results which are specific to the TASEP and 
give explicit formulae for the generalized hamiltonians Hk- The calculations 



Golinelli, Mallick — Commuting operators for TASEP 



8 



leading to these expressions are carried out in detail in Appendix [Bi 

The t(0) operator appearing in equation ffTBT) is the translation operator 
on the ring and is defined as 

*(0)|ri,r 2 ,...,r £ ) = |r 2 , . . . , t l , n) . (18) 

The operator H\ given by 

L 

H 1 = t'(0)/t(0) = Y,M i = M, (19) 
i=i 

is precisely the Markov matrix M which thus belongs to the commuting 
family generated by t(X). 

All the iffc's can be explicitly calculated. By using the "ring-ordered 
product" defined in section 1231 we find in Appendix IB1 that for 1 < k < L, 

H k = Yl 0(M h M i2 ...M ik ) . (20) 

l<il<i2<---<ik<L 

In particular, we have 

H 2 = J2 OiMiMj) , (21) 
i<*<j<£ 

and according to equation (TT%|) . 

H L = |1,1,...,1)(1,1,...,1| . (22) 

For k < L, all the terms in Hk are products of k jump matrices which, because 
of ring-ordering, correspond to k different particles jumping simultaneously 
one step forward. Thus, has a non- vanishing action only on configurations 
with at least k particles. 

In the case with L = 4, for example, the generalized hamiltonians are 
given by 

H x = .U, • M, ■ .U ;i • .U, M . (23) 

//•_, = .U|.U, + M 2 M 3 + M3M4 + .\/;.\/| + + .Uj.U; , (24) 

//.; = M1M2M3 + M2M3M4 + M3M4M1 + M4M1M2 . (25) 

Using Eqs. ( TTfiT) and (j2T?l) . we conclude that the generating function of the 
Hk is given by 

t 9 (X) = H = 1 + g A fc /f fc = f + AM 4 ) j . (26) 
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Although the operator Hi is the Markov matrix M of the TASEP, we empha- 
size that when k > 2, H k cannot be interpreted as a Markov matrix because 
it contains negative non-diagonal matrix elements. However, we shall now 
prove that the matrix t(X) is the Markov matrix of a discrete time process 
when < A < I. 

3.2 Action of the transfer matrix on a given configu- 
ration 

We describe now the action of t(X) on a given configuration |r 1; r 2 , . . . , r L ). 
Using equation ( 1261) . we observe that 



t g (X)\T h T 2 ,...,r L ) = 



( L \ 

Hil + XMi) 



\ti,t 2 ,...,t l ), (27) 



where the product runs only over occupied sites. This expression shows that 
the action of t g (X) is factorized block by block. We consider first the simple 
block |0P) (the notation |01 p ) means that a hole is followed by p particles). 
We have 

v 

t(A)|01 p > = t,(A)|P0> = (1 + XMi) . . . (1 + AM p )|P0> = J2fkJ± k W k ) , 

with f 0tP = X p and f k>p = (1 - X)X p ~ k for 1 < k < p. (28) 

More generally, for a configuration C of the form |0' 31 1 P1 92 1 P2 . . . gs l Pa ) with 
Pi, Qi > 1) we obtain 

t(x)c = t g (X)\o 9i ~ 1 i pi . . . o q n Ps o) = 

lo^- 1 )® ( f kl:Pl \l kl 01 pl ~ kl ) ) ® lo" 2 - 1 ) ® . . . 

\A:i=0 / 
' Ps \ 

X)/*.*.|l*"01*--*-> ) . (29) 



Except for the full configuration |1 L ) for which t(A)|l L ) = t g (X)\l L ) = 
(1+A L ) |1 L ) and the void configuration |0 L ) for which t(X) |0 L ) = t g (X) \0 L ) = 
|0 L ), any configuration C has at least one particle and one hole. By using 
the translation operator t(0) that commutes with t(X), it is possible to bring 
C to a form to which equation (|29il can be applied. 



Golinelli, Mallick — Commuting operators for TASEP 



10 



We notice that t(l) is the identity operator. Consequently t 9 (l) is the 
forward translation operator, t g (l) = t(0) _1 . 

We illustrate these results with an example of 3 particles on a ring of 5 
sites: 

t(A) |10101) = t g (X) |01011) = 

(1 - A) 2 |01011) + A(l - A) 100111) + A(l - A) 2 |11010) 

+A 2 (1 - A) |10110) + A 2 (l - A) |11001) + A 3 |10101) . (30) 

By considering the action of the operators t(X) and t g (X), we remark that 
for < A < 1, fk ;P > and that J2k fk,p = 1- The quantities fk ;P can thus 
interpreted as probabilities. The operators t(X) and t g (X) are then Markov 
matrices of discrete time exclusion processes with parallel dynamics, in which 
different holes can jump simultaneously through clusters of particles. 

With t(X), a hole located on the left of a cluster of p particles can jump a 
distance k in the forward direction, 1 < k < p, with probability X p ~ k (l — A). 
The probability that this hole does not jump at all is A p . 

With t g (X), a hole located on the right of a cluster of p particles can 
jump a distance k in the backward direction, with probability A fc (l — A) for 
1 < k < p, and with probability A p for k = p. The probability that this hole 
does not jump at all is 1 — A. 

The t g (X) Markov process is equivalent to a 3-D anisotropic percolation 
model and a 2-D five- vertex model (Rajesh and Dhar 1998). It is also an 
adaptation on a periodic lattice of the ASEP with a backward-ordered se- 
quential update (Rajewsky et al. 1996, Brankov et al. 2004), and equiv- 
alently of an asymmetric fragmentation process (Rakos and Schiitz 2005). 
Consequently Markov matrices of these models on a periodic lattice form a 
commutating family. 

3.3 Invariance properties of the transfer matrix 

We describe here the symmetries of the transfer matrix t(X) and of the op- 
erators iffc. Translation invariance is obvious because t(0) is the translation 
operator and commutes with t(X) and H^. From Eqs. (T2T)} |2"B"1) . we observe 
that t(X) and Hk conserve the number n of particles because each jump 
matrix M, does so. For a given value of n, t(X) is a polynomial of degree n. 

The Markov matrix M is symmetric under lattice reflection R (obtained 
by exchanging sites i and L — i + 1) followed by particle-hole conjugation C 
(Golinelli and Mallick 2004). This CR symmetry acts on a configuration as 
follows 

CR\t 1 ,t 2 , . . . ,t l ) = |1 -tl, ■ ■ ■ , 1 -r 2 , 1 -ri) . (31) 
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The CR symmetry does not commute with the translation operator t(0) 
because CR t(0) = t(O)" 1 CR. The following property 

CR M { CR = M L _i , (32) 

implies that CR is a symmetry of the Markov matrix, i.e., (CR)M(CR) = 
M. However, CR is not a symmetry of Hk for k > 2 because the orientation 
of matrices along the ring is inverted by Eq. ( |32i) . More precisely, Hk and 
t(X) are transformed as follows 

H k = CR H k CR , t(X) = CR t(X) CR , (33) 

where H k and t(X) are given by formulae similar to Eqs. ( 1201 1261) but with an 
antz-ring-ordered product O instead of the ring-ordered product O. With the 
H k operator, k different holes jump simultaneously one step backward. For 
k > 2, one can verify that the action of Hk and of Hk on a given configuration 
are different. 

The CR symmetry allows us to construct two different families t(X) and 
t(X) of commuting operators, i.e., [t(A),t(z/)] = and [t(A),t(z/)] = 0. Both 
families contain the Markov matrix M = H\ = Hi. However, t(X) and t(i>) 
do not commute with each other for generic values of A and v. 



4 Connected Operators 

In the previous section, we have defined a set of commuting operators, the 
generalized hamiltonians H k , that act on k different particles. However, 
these actions are generally not local because they involve particles with ar- 
bitrary distances between them. Moreover, as can be seen from Eq. ( 1201 ). 
the number of terms in H k for a large system of size L grows as L k /k\. In 
statistical physics, quantities that are local and extensive are preferred. Such 
"connected" (or local) operators are usually built from the logarithm of the 
generating function (Liischer, 1976). Therefore, for k > 1, we define the 
connected hamiltonians F k as follows : 

oo ^ k 

ln^A)^^. (34) 

k=l 

The Ffc's can be expressed from the H k s using definition ( |26l) . Be- 
cause the H^s are commuting matrices, the F^s are also a set of commut- 
ing operators and moreover they commute with t(A) and with all the H^s, 
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i.e. [F k ,Hj] = 0. Consequently the F k s can be calculated with the usual 
moment s-cumulant s t r ansf or mat ion , 

fc-i 

F k = kH k -J2 F i H k-i (35) 
i=i 

which is obtained from the derivative of lnt ff (A). 

We now show that lni s (A) and the F k s are linear combinations of con- 
nected words, i.e., words which cannot be factorized in two (or more) com- 
muting words. Consider a word W of lni 9 (A) made of jump matrices Mj 
with % G X C {1, 2, . . . , L}. This word must also appear in lntj(A) with 

t I (\) = o(^l(l + \M i ) ) j . (36) 

Assume that the set of indices X can be split into two disjoint subsets, Xi and 
X2, such that [M a , Mb] = for all a G Xi and all b G X2. Then the ring-ordered 
product in equation f|36|) can be factorized in two non-connected products 
and we have 

lntx(A) =lnt Zl (A)+lnt Z2 (A) . (37) 

Therefore W must be made of jump matrices with indices % all belonging 
either to X\ or to X 2 . Applying this reasoning recursively, we deduce that W 
is connected. We emphasize that connected words remain connected after 
the use of the simplification rules §7S§§ . 

4.1 Calculation of Fk for small k 

We first remark that Eq. (134p defines an infinity of operators F k but we have 
seen that there are only L operators H k for a system of size L. Therefore, the 
Xfc's are not all independent and the knowledge of the Fi, . . . , Fl is formally 
sufficient to generate all the F k . Consequently, when we consider F k in the 
following formulae, we assume implicitly that the system is sufficiently large 
to have k < L. The operator F k is the fcth order term in the expansion of 
\ntg(X), given by equation (l33|l . After using the relations and ([B]), F k is 
found to be a linear combination of words of length j, with j < k. 
For k = 1, Fx is the Markov matrix M, 



L 

F l = Hi = M = Mi 



(38) 
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For k = 2, we have 

L 

F 2 = 2H 2 - Hi = Y^iWu M i+ i] + Mi) , (39) 
i=i 

where we use the convention M^l = Mi due to periodic boundary conditions. 
The operator F 2 is indeed connected : all non-connected terms in 2H 2 — H\ 
of the type MjMj with \% — j\ > 2 cancel one another and there remains 
only words of the type MjMj +1 and M i+1 M i: involving the adjacent bonds 
1) and (i + l,i + 2). 
After an explicit calculation, we find the following formulae for F 3 , F4 
and F 5 : 

F 3 = 3H 3 — 3H 2 H 1 + Hf 

L 

= ^([[Mj, Mi +1], M i+2 ] + MiM i+ i — 2M i+ iMi + Mi) ; (40) 
i=i 

F 4 = AH 4 - 4# 3 # 1 - 2Hl + 4H 2 Hf - H\ 

L 

= ^2 { [[[Mi, Mi +1], M i+2 ], M i+3 ] 

i=l 

+MiM l+1 M t+2 - 2{M i+1 MiM i+2 + M i+2 MiM i+1 ) + 3M i+2 M i+1 Mi 
+MiM i+1 - 3 .\ /,.,.!/, + M } ; (41) 
F 5 = 5# 5 - 5# 4 #i - 5# 3 #2 + 5F 3 #i + 5# 2 2 #i - 5ff 2 #i + Hi 

L 

= ]T{[[[[M,M +1] ? M i+2 ] , M i+3 ] , M i+4 ] 

i=l 

+M i M i+1 M i+2 M i+3 

-2(M i+1 MiM l+2 M i+3 + .\/,.,.\/,,\ /,.;.!/,., + M l+3 MiM i+1 M i+2 ) 
+3(M i+2 M i+1 M i M i+3 + M i+3 M i+l MiM i+2 + .U, . :i .1/, . .U,,U, . : ) 
-4Af <+3 M H . 2 Af i+ iM i 

+MiM t+l M t+2 - 3(M i+1 MiM i+2 + M i+2 MiM i+l ) + (i.U, .•>.!/, . : U 
+MM+X - IM.iM + M } . (42) 

As expected, is made only of connected words. We notice the following 
remarkable property from the expressions (I40H42I) : the words of length j in 
Ffc are always a permutation of j consecutive matrices, M i? M i+1 , . . . , Mj +J -_i, 
without repetition. For example, the expression (j4*Tl) of F 4 does not contain 
the word Mi + iMjMj+2Mj + i. This property of has been verified explicitly 
for k < 10. 
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4.2 A formula for the connected operators 

We have written a computer program that gives the expressions of the F^s 
for small values of k (up to k = 10). This leads us to conjecture a general 
formula for F k valid for arbitrary k. In order to write this general formula 
we need to define some notations. 

4.2.1 Simple words 

A simple word of length j is defined as a word M ff mM ff (2) . . . M a u\, where 
cr is a permutation on the set {1, 2, . . . , j}. For example, there is a unique 
simple word of length 1, noted W\ = Mi and two simple words of length 2, 
W 2 {1) = M X M 2 and W 2 {0) = M 2 M X . For j > 2, the commutation rule flSJ im- 
plies that only the relative position of Mj with respect to Mi±\ matters : the 
number of simple words of length j is therefore much smaller than j ! . In fact , 
any simple word Wj is uniquely characterized by (s 2 , s 3 , ... , Sj) where Sj = 1 
if Mi is written to the right of Mj_i in Wj and Sj = otherwise. Therefore, 
there are 2 J_1 simple words of length j and we note them Wj(s 2 , S3, ... , s^). 
Simple words obey the recursive rule: 



The set of simple words of length j will be called Wj. 

For a simple word Wj, we define u(Wj) to be the number of inversions in 
Wj, i.e., the number of times that Mj is on the left of Mj_i : 



By definition, < m(W}) < j — 1. For example, we have ^(l, 0,1,0) = 
M 5 M 3 Mi M 2 M 4 and m(W 5 (1, 0, 1, 0)) = 2. 

Using these definitions, the "nested" commutator that appears in the 
expressions ( 13911421) can be written for general k as : 



W j (s 2 ,s 3 ,...,s j ^ 1 ,l) 
W j (s 2 ,s 3 ,...,s j ^ 1 ,0) 



V^-i(s 2 ,s 3) • • -,Sj-i) Mj 
MjWj-i(s 2 , s 3 , ... , Sj-i) . 



(43) 
(44) 




(45) 



[[...[[M 1 ,M 2 ],M 3 ],...},M k ]= (-lT m W, 



(46) 



wew k 



where ^2 



is equivalent to writing J2 S =0 • • • S 



1 



wew k 
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4.2.2 Conjectured general formula for F k 

We have calculated the exact expressions of the connected operators up to 
Fiq and we have noticed that in F k all simple words W of length j < k 
appear with the sign (— and with a coefficient given by the binomial 
coefficient (* ~°k^ ) ■ Therefore, for k < L, we conjecture the following 
general formula for F k : 

^ = T E ^(-ir^f k - j ^ w) ) w , (47) 

j=l W&Wj ^ 3 ' 

where T is the translation-symmetrizator that acts on any operator A as 
follows : 

L-l 

TA = Y,t{0) 1 At^O)-' . (48) 

The presence of T in equation (|4T|) insures that F k is invariant by translation 
on the periodic system of size L. 

We also verified that for j + k < 11 the conjecture (HTj) gives [F k , F 3 ] = 0. 
We emphasize that because of the special expression (1221) of Hl the expression 
(I4T1) of F k is valid only for systems with length L > k. 



4.3 Action of Fk on a configuration 

In this section we describe the action of Fk, as given by the formula (1471) . on 
an arbitrary configuration C = |ti, t%, . . . , Tl). We first define an operator 
A, that we shall call the 'Antisymmetrizator', by describing its action on a 
configuration. The antisymmetrizator A acts on a bond as follows : 

A|01) = |01> — |10> , (49) 
and A \Ti,T i+1 ) = {r^n+i) , for r» ^ and r i+l ^ 1 . (50) 

More generally, the action of A is given by : 

A |lP 1 0' ll+1 l P2 + 1 ?i2 + 1 l P2 + 1 Ofrr-l+l]Pr+lQ?V\ — 

|F 1 hl ) ® A |01) ® (1^0^) (8) A |01) ® . . . ® A |01) ® IF-O^) , (51) 
where h iy pi > 0. 

Consider now a simple word Wj(s2, s 3 , . . . , Sj) acting on a system of size 
L > j. This operator affects only the sites 1, 2 . . . , j, the sites j+2, . . . , L 
being spectators. We show in Appendix O that 

Wj(s 2 , s 3 , . ■ ■ , Sj) |n, r 2 , . . • , tl) (52) 
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if and only if 

n = 1, r j+1 = and r 2 = s 2 , r 3 = s 3 , . . . , = s 3 - . (53) 
If this condition is satisfied, the action of the simple word is given by 

Wj(s 2 , S 3 , . . . , Sj)\l, S 2 , . . . , Sj, 0, T j+2 , ...,T L ) = 

A |0, s 2 , . . . , Sj, 1) <g> |r i+2 ,...,T L ) , (54) 

where A is defined in equation ( l5Ti) . Thus, a word acts only on specific 
configurations. From this remark, we can derive a formula for the action of 
Ffc on a configuration C. From equation fT4T|) . we first observe that only one 
specific word G VVj has a non zero action on a given configuration C : 



E^r^C^V 11 '" 2 '---'^' ' 



■ ■ ■ ,Tl) 



(-l) u r^)A\0,T 2 ,...,T j7 l) ® |r i+a ,...,7i), (55) 

where m = J^^ =2 (l — r ?) is the number of holes in C between sites 1 et j '• + 1. 
We emphasize that u is a function of C only. Now, according to equation (jUj), 
we have to take a sum over j and apply the translation-symmetrizator T . 
This amounts to considering all possible jumps from an occupied site i to an 
empty site m with j = m — i < k. We thus obtain 

F k \n,r 2 , . . . ,r L ) = 

L i+k 

-. . „,; / k+i-ii) + iHi.i)i) \ ^ 



1^ / \ 



i=l m=i+l 

Tj=l T m =0 



In, . . . ,7i_i) g) A|0,r i+1 , . . . ,r m _i, 1) ® |r m+ i, . . . ,r L ) , (56) 



where u(i,m) = YlT=i+i^ ~~ T r) * s ^ ne number of holes in C between sites i 
and m (we recall that sites are defined modulo L). 

The action of can be described as follows. Each particle, starting from 
an occupied site, can make all possible jumps of length j < k to a vacant 
site. Each jump has a sign and a weight : the sign is given by (—1)" where 
u is the number of holes overtaken by the particle between its initial and its 
final position; the weight is a binomial coefficient that depends only on u and 
k-j. 
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5 Conclusion 

The algebraic Bethe Ansatz technique allows to construct a family of opera- 
tors that commute with a given integrable hamiltonian. For the totally asym- 
metric exclusion process, this procedure has enabled us to define a family of 
generalized operators, local and non-local, that commute with the Markov 
matrix. The properties of these operators have been derived by using the 
TASEP algebra (JTHHJ) and their actions on the configuration space has been 
explicitly described. In particular, we have found a combinatorial formula for 
the connected operators valid at all orders. This formula has been verified 
for systems of small size but a mathematical proof remains to be established. 

It would be of interest to extend our results to the exclusion process with 
forward and backward hopping rates. Because the symmetric exclusion pro- 
cess is equivalent to the Heisenberg spin chain, the generalized hamiltonians 
would correspond to integrable models with long range interactions. Explicit 
formulae for the connected conserved operators associated with the Heisen- 
berg spin chain are known only for the lowest orders (Fabricius et al., 1990); 
no general expressions for these spin chain operators have yet been found. 
We believe that the expression given in the present work, equation (H7j), that 
is valid at all orders, may shed some light on this issue. 

Finally, we hope that the family of commuting operators studied in the 
present work will help to explain the spectral degeneracies found in the 
Markov matrix and to unveil hidden algebraic symmetries of the exclusion 
process. 
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A Construction of the transfer matrix 

In this appendix, we use the algebraic Bethe ansatz method to construct the 
transfer matrix of the TASEP, a family t(X) of commuting operators acting 
on the configuration space TLl- 

Generalized jump operators 

In section [221 we have defined Mi, the jump operator from site i to site i + 
More generally, for two different sites a and b, we define P a ^i the permutation 
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operator between sites a and b by 

P ajb \ . . .r a , . . . r 6 , . . .} = | . . .t 6 , . . . r OJ • • •} , (57) 
and M a> b, the jump operator from a to b, by 

M a , b = (P a>b - 1) a«aj 0) (58) 

where a- r ^ = |rj)(rj| is the projector on the subspace of configurations with 
site % in state r. The operators M a ^ and P aj b act non trivially only on the 
subspace V a <S> V b and are the identity operator on all spaces Vi for i different 



from a and b. The relations ([7H9]) now become 

Ml b = -M a , b (59) 

M a , b M hiC M af> = M btC M a , b M bjC = (60) 

[M a , b ,M c J = (61) 



where a, 6, c and d are different sites. Equation (1581 allows to define a totally 
asymmetric exclusion process on an arbitrary graph, with one jump matrix 
M atb for each directed edge (a, b) of the graph. Consequently Mi is just a 
simplified notation for M^+i when the graph is a ring. 

As the main problem is the non-commutativity of the neighboring jump 
operators Mj and M i+1 , the key step consists of finding operators which 
verify a quasi-commutation rule, the Yang-Baxter equation. Such operators 
are given by 

C a>b (X) = P a , b (l + \M a;b ), (62) 

where a and b are two given sites, and A is a number (the spectral parameter). 
The C a fi satisfy the Yang-Baxter equation (for a derivation see, e.g., Golinelli 
and Mallick 2006b) : 

C a>b {v)£ c>b {\)£ C;a {ij) = C C! a(n)C c>b {\)C a>b (v) if A = p + v - [iv. (63) 
The monodromy matrix T(X) 

To the L physical sites (z = 1, . . . , L), we add an auxiliary site with label 0. 
The extended configurations are noted as |tq) <g> \ti, . . . , tl), with G {0, 1} 
for i — 0, . . . , L, and the extended 2 L+1 dimensional state space is given by 
Vo <8> H-l- I n order to distinguish the spaces on which operators act, we note 
with a "hat" ( " ) the operators acting on the extended space Vq ® 'Hl, an d 
without hat those acting on the physical space Hl- 
We define the monodromy matrix T(A) by 



f (A) = A,o(A) £ 2i0 (A) . . . £ ij0 (A) 



(64) 
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The matrix T(A) acts on the extended space V (g> TCl- We now consider 
two auxiliary sites et 0' and two monodromy matrices T (A) and T >(u) 
acting on the space V$ ® Vqi ®TCl- Using equation ( l63l) and the fact that 
\Ci<a(X),Cj<oi(ii)] = for i ^ j, we deduce that T (A) and T /(/i) also satisfy 
the Yang-Baxter relation : 

Co>,q{v) To (A) To/(/i) = T /(/x) T (A) £ ',o(^) if X = + u - /jlu . (65) 

Using definitions ([6"2"j |6"^|) we find for A = 0, 

f (0) = A,o A,o • • • Pl,o. (66) 

The explicit action of T(0) on an extended configuration is then 

r(0)(|ro>®|ri,r 2 ...,r L )) = |n> ® |r 2 . . . , r L , r ). (67) 

It turns out that T(0) is the translation operator which causes a left circular 
shift of the sites, including the auxiliary site 0. 

In Eq. ( I64p for a generic A, we can "push" the permutation operators P^o 
to the left using the relation Mj-^o A,o = A,o and obtain 

f (A) = f (0)(1 + AM 1)2 )(1 + AM 2i3 ) . . . (1 + AM L _ ljL )(l + AM L , ) . (68) 

The operator T(A) is a polynomial of degree L, 

f(A)=f(0) ^l + X> fc T^) ( 69 ) 

where the T\'s, that act on Vq ® Hl, are given by 

f fc = M iuil+1 M i2 , l2+1 . . . M lk , ik+1 (70) 

l<ji<«2<...<ife<I/ 

for 1 < k < L, with the convention Ml,l+i = ^l,o when = L. Hence, 
the operator T^ represents the simultaneous jumps of k different particles 
initially located on the physical sites. In particular, T\ is the Markov matrix 
of the TASEP on the open segment (1, 2, ... , L, 0). 

The trace over the auxiliary space 

As the operators defined above act on the extended space Vq ® TCl, we will 
use the partial trace tro over the auxiliary space Vq to obtain operators acting 
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only on the physical space TCl- Any operator A acting on V <S> TCl can be 
uniquely written as 

1 

i= £ (WW)® (7i) 

where A(to, t' q ) is an operator acting on TCl- The partial trace is defined as 

i 

tr i=E A (ro, r ) = A(0, 0) + A(l, 1) , (72) 

r =0 

and the action of tr A is given by 

i 

tr A|C> = 5~>o|i(|ro>®|0), (73) 

r =0 
1 

(C'|tr i=E(< r o|®(C / |)i|To>. (74) 

r =0 

Another property of the trace that we shall need is the following. Consider 
an operator X that acts only on TCl', this operator can thus be written as 
X = 1 ® X. Then for any A acting on Vq <E> Tii we have : 

tr (iX) = tr (i)X, (75) 
tro(li) = Xtr (i). (76) 

The transfer matrix t(\) 

The transfer matrix t(\), which acts on the physical configuration space H,l 
is defined by 

t{\) = tr T(A) . (77) 

The operators £j(A) and the monodromy matrix T(A) conserve the number 
of particles in the extended space (physical space plus the auxiliary site). As 
the auxiliary trace tr operation keeps constant the number r of particles 
on the auxiliary site, it keeps the number of particles in the physical space 
constant too. Hence by construction, the transfer matrix t(X) conserves the 
number of particles. 

We now multiply the relation (!B3|) by C ^(v) on the left and take its 
trace tr ,o' over the two auxiliary sites and 0'. Because Cq} q {v) acts only 
on and 0', we can use that tr ,o' is cyclic with respect to C^ Q (i>) and thus 

tro.cPoCA) Tv{p)\ = tr , 'M)u) T (A)] . (78) 
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Using ([77]) and the relation tr ' = tr tr /, we obtain 

t(\) t(n) = t(\) t(ji) . (79) 

The Yang-Baxter equation (["""*]) thus implies the commutativity of the trans- 
fer matrices. 

B Calculation of the hamiltonians 

We derive here the expression ([2*01) for H k . Following Eqs. ([TOj 1001 177]) . we 

obtain, for 1 < k < L, 

t(0) = tr [f (0)] and H k = i(O)- 1 tr [T(0)T fc ] • (80) 

We shall now perform the trace tr over the auxiliary space. 

We first calculate t(0) : for a given configuration |n, T2, . . . , n) of the 
physical sites, we obtain using Eqs ( [731) and ( [771) 

i 

£(0)|ti,t 2 , • . .,T L ) = ^ (r |T(0)( |r ) ® |n, r 2 , . . . , r L » . (81) 

T =0 

As T(0) is the translation operator on the extended space, we obtain 

i 

t{o)\n, t 2 ,...,t l ) = J2 ( r o\ ( In) ® |r 2 , . . . , n, r )). (82) 

r =0 

In the auxiliary space Vo, we have (nln) — ^r ,ri and then 

*(o)|Ti,r 2 ,...,Tfc) = |r 2 , . . . ,r L ,n) . (83) 

Therefore t(0) is the translation operator on the configuration space. 

We now evaluate H k for 1 < k < L — 1. According to Eqs. ([701 150]) . any 
term that appears in T k is made of k jump operators with k < L. Thus, 
such a term W can always be written as DF with 

D = M iuil+ i Mi 2 ,i 2 +i ■ ■ ■ M iriir+1 with l<i 1 <i 2 <...<i r <u-l, 
F = M ir+uir+1+1 . . . M ik ,i k+ i with u + 1 < i r+1 < ... <i k < L. (84) 

The index u is such that the matrix M U)U +i does not appear in W . Therefore, 
all the traces that we have to calculate are of the type tr (T(0)DF) with 
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[D, F] — 0. Besides we notice that D acts only on Ti^. Therefore, we have, 
using equation flTol) : 

tr (r(0)DF) = tv (f(0)FD) = tr (f(0)F)D , (85) 

with .D = Mi x Mi 2 . . .M ir . The operator F can not be extracted from the 
trace because it acts on the auxiliary site if if. = L in equation ( j84l ). 
However, recalling that T(0) is the translation operator on the total space 
Vo <8> TCl, we can write 

F (0)F = F' f (0) with F' = M v+1 ^i, v+1 • • • M ik - 1)ih . (86) 

The fact that u > 1 and ik < L ensures that F' acts only on and as such 
can be written as 

F' = 1 ® F' with F' = M ir+1 _x . . . M ifc _x . (87) 

We now use the property ( 1761) and write equation ( 1851) as follows 

tr (f (0)/)F) = tr (F'f (0))D = F'tr o (f(0))D = F't(0)D . (88) 

Using the fact that t(0) is the translation operator on the configuration space, 
we write 

F't(0) = t(0)F with F = M ir+1 ...M ik , (89) 

and conclude that 

tr (f (0)DF) = t(0)FD = t(0)O (DF) , (90) 

where O () is defined in section [2~51 This proves the general formula (1201) . 

To be complete we need to calculate the operator of the highest degree 
Hl- The operator 

f L = M 1>2 . . . M L _ 1>L M Lfi (91) 

involves jumps from all physical sites : it can not be split as described in 
equation (I51I) . We have T L \C) = for all configurations C, unless C = 
|0) ® |1, 1, • • • 1). After a short calculation, Eq. (ISO]) leads to that 

# L = |1,1,...,1><1,1,...,1| , (92) 

which is the projector on the "full" configuration (all sites are occupied) in 
agreement with equations (tHJ) and (1221) . 
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C Derivation of Eq. ([52HMj) 



In this appendix, we prove equations (15241541) by induction on the size j of 
the simple word W = Wj(s 2 , S3, ... , Sj). We shall simplify the notations by 
writing the action of W on the sites 1, 2 . . . , j, j + 1 (the sites j + 2, . . . , L 
being spectators). 

For j = 1, the only word is W\ = M\ and equations (15211541) are satisfied : 



M 1 \t 1 ,t 2 )^0 iff n = l, r 2 = (93) 
and Mi |1,0) = |0, 1) - |1,0) = A |0,1). (94) 

For j > 2, we shall calculate the action of the word W on the configuration 
C = \ti,t 2 . . . Tj+i). We must distinguish two cases Sj = 1 or 0. 



Case Sj = 1 

We can write W = W Mj where W = Wj-±(s 2 , . . . , Sj_i) is a simple word 
of length j — 1 and we have 

w |n . . . n-i, n, n+i) = W M,- |n • • • T jt r j+1 ) (95) 

This action vanishes unless Tj = 1 = Sj and n +1 = 0. In that case we have 

win... Tj . x , 1, o) = W |n ■ ■ ■ n-i> o, l) - W |n ■ ■ ■ n-i> i> °) • ( 96 ) 

We can now use the induction hypothesis : the second term on the r.h.s. 
always vanishes (because Tj = 1); the first term on the r.h.s. does not vanish 
if n — 1 an d r 2 = s 2, ■ ■ ■ , Tj-i = Sj^i- Therefore, the action of W on C does 
not vanish if and only if C = |1, s% . . . 1, 0) and is given by 

W\l,s 2 ... 1, 0) = W |1, s 2 . . . Si _i, 0,1) = A |0, 8 2 ... sj-!, 1, 1) , (97) 

where we have used the induction hypothesis to evaluate the action of W 
(we recall the site number j + 1 is spectator for W). Equations (152H54P are 
thus proved for the case Sj = 1. 



Case Sj = 

We now have W = Mj W where W is defined as above. Therefore 

W n_ 1; Tj, Tj+1 ) = Mj W \n--- n_i, r 3 , Tj+1 ) . (98) 

The induction hypothesis implies that the action of W does not vanish if 
and only if n — 1, r 2 — s 2 , ■ ■ ■ , n~i = s^-i, Tj = = Sj, the site (j + 1) 
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being spectator for W. Besides, the action of Mj on the bond + 1) is 
non-trivial only if r J+1 = 0. Therefore, we have C = |1, s 2 • • • Sj-i, 0, 0) and 

W |1, s 2 . . . 0, 0) = Mj (A |0, s 2 ... 1) <g> |0) ) . (99) 

The action of A on the r.h.s. of this equation is given by 

A |0, s 2 • • ■ Sj-2, 1, 1) = A |0,s 2 ... 3^2,1) ®|1> if 3^1 = 1, (100) 
A |0, s 2 • • ■ Sj_ 2 , 0, 1) = A |0, s 2 . . . Si _ 2 ) <8> A|0, 1) if Sj -_i = . (101) 

Thus, we obtain, if s_j_i = 1, 

W |l,s 2 -..Sj-2, 1,0,0) = A |0, s 2 . . . Sj_ 2 , 1) ® Mj\l, 0) 

= A |0, s 2 . . . Sj_ 2 , 1) ® A|0, 1) 
= A |0, s 2 . ..Sj-2, 1, 0, 1} , (102) 

and if = 0, 

W |l,s 2 ...s 3 -_ 2 , 0,0,0) = A |0,s 2 ...s 3 -_ 2 ) ® M,-(A|0,1) ® |0)) 

= A |0,s 2 . . .s 3 -_ 2 ) ® |0) ® A|0,1) 
= A |0,s 2 ... Sj-2,0, 0,1), (103) 

which completes the proof of Eq. ( 1541) . 
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